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Abstract 

The physics of fluid interfaces between domains of different magne- 
tization in the ordered phase of the 3D three-state Potts model is stud- 
ied by means of a Monte Carlo simulation. It is shown that finite-size 
effects in the interface free energy are well described by the capillary 
wave model at two loop order, supporting the idea of the universality 
of this description of fluid interfaces in 3D statistical models. 
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1 Introduction 



It is well known that 3D spin systems on finite volumes show domain walls 
separating coexisting phases which behave as fluid interfaces between the 
critical and the roughening temperature. The finite-size effects in the free 
energy of a fluid interface are dominated by long-wavelength fluctuations 
and a correct physical description of the critical properties of the surface 
cannot neglect their contributions Q. 

While below the roughening temperature, where the interfaces are almost 
rigid, a microscopical approach can be taken (see Q and references therein), 
above it one is forced to assume an effective model describing the collective 
degrees of freedom of rough interfaces. 

The capillary wave model (CWM) || in its simplest formulation, which 
we will follow, assumes an effective hamiltonian proportional to the area of 
the surface. 

It has been recently shown [Q] that rather strong finite-size effects, de- 
pending on the shape of the lattice, are well described in terms of the one- 
loop or gaussian approximation to the CWM: its predictions have been tested 
with high accuracy in the scaling region of the 3D Ising model. These one- 
loop corrections depend only on one adimensional parameter, namely on the 
asymmetry z = R/T of the transverse sizes of an elongated lattice RxT x L 
(L » R, T) with periodic boundary conditions taken in all directions |], ||, [|] . 

As it has been already pointed out ||], higher order corrections to the 
gaussian model can be taken into account to verify the CWM beyond the 
one-loop approximation. These higher order corrections can provide a more 
stringent test of the CWM model: in fact, while many different effective 
hamiltonians reduce to the gaussian form at one-loop level |?J ; they differ 
in the form of two- and higher-loop corrections. In this paper we apply the 
CWM in the two-loop approximation to the finite-size behavior of order- 
order rough interfaces in the 3D three-state Potts model. 

The two-loop contributions do not depend only on the asymmetry pa- 
rameter z but also on an adimensional parameter proportional to the min- 
imal area of the surface, namely o~RT, a being the (reduced) interface ten- 
sion. An important consequence of this fact is that the well known finite-size 
behavior ||] of the energy splitting E occurring between vacua on finite vol- 
umes, has no longer, at two-loop, the classical functional form 

E d (R) oc e-° R2 (1) 

not even for symmetric (R = T) lattices. This must be contrasted with what 
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happens at one-loop level, where, due to the scale-invariance of the one- 
loop contributions, one has to consider asymmetric lattices to find deviations 
from the classical functional form Q. 

To verify the CWM we study the finite-size behavior of rough order- 
order interfaces of the 3D three-state Potts model by means of Monte Carlo 
numerical simulations. 

This model is of great interest because of its well known connection with 
4D SU(3) pure Yang-Mills theory at finite temperature Q. In fact one can 
assume an effective action of a 3D spin model with short-range, ferromag- 
netic interaction to describe the finite temperature deconfinement transi- 
tion of QCD in the limit of infinite quark masses |flO|| . For this reason the 
Potts model has been extensively investigated and, in particular, the prop- 
erties of order-order and order-disorder interfaces have been already studied 
[11, 12, 13], one of the main goals being the evaluation of the corresponding 
interface tensions. 

In this paper we show that the knowledge of the functional form of the 
fluid interface free energy on finite volumes, including the capillary wave 
contributions, enables one to estimate the order-order interface tension with 
high precision. Moreover, the good agreement we find with the theoretical 
predictions of the CWM model at two-loop order strongly supports the idea 
of the universality of this description of rough interfaces in 3D statistical 
models. 

The paper is structured as follows: in Sec. 2. we describe the CWM 
and its one-loop and two-loop approximations. In Sec. 3. we compare the 
results of MC simulations with the CWM predictions. Sec. 4 is devoted to 
some concluding remarks. 



2 The capillary wave model 

According to the CWM, the interface between two domains of different mag- 
netization in the ordered phase of a 3D spin system, above the roughening 
temperature, is described by the partition function 

Zaw = J [Dx]exp{-aA[x]} , (2) 

where the single- valued function x(r, t) describes the displacement from the 
equilibrium position of the interface, a is the reduced (order-order) interface 
tension and A[x] is the area of the interface 
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It should be mentioned that (|3|) coincides with the Nambu-Goto string action 
in D = 3 in a particular gauge. Eq. (||) is not expected to be the exact action 
describing fluid interfaces but at least the dominant contribution^]: as we will 
show this is indeed the case. 

To compare the predictions of the CWM with numerical results from 
Monte Carlo (MC) simulations, we have chosen 3D lattices of R x T x L 
sites, with L S> R,T, and periodic boundary conditions in each direction. 
This particular choice of the lattice shape allows one to consider only in- 
terfaces orthogonal to the elongated (L) direction, the probability of having 
interfaces orthogonal to the other directions being negligible. The 2D field 
x(r, t) is therefore defined on the rectangle (r, t) 6 [0, R] x [0, T] with opposite 
edges identified, i.e. on a torus. 

The partition function (||) cannot be computed exactly, but it is possi- 
ble to express it as an expansion in powers of the adimensional parameter 
{(tRT)^ 1 : the two-loop expansion of Z cw can then be written as 



Z CW (R,T) oc e 



jRT 



y(l loop) 



■Zf lo °^ (R,T) 



(4) 



The one-loop contribution (namely the gaussian approximation), obtained 
retaining only the quadratic term in the expansion of (^), is nothing else 
than the exact partition function of a 2D conformal invariant free boson on 



a torus of modular parameter r 
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while the two-loop term can be calculated perturbatively expanding (J3T) at 



the next-to-leading order [15] 

1 
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"'The problems arising in the quantization of non-critical strings are known to disappear 
asymptotically at large distances [Jwj. 
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The two functions r] and E2 appearing above are respectively the Dedekind 
function and the second Eisenstein series: 



00 



r,( T ) = q 1/24 Y[(l-q n ), Q = exp(2vrir) 
n=l 

00 n 

E 2 (t) = l-2AY^— . 

^ 1-Q n 



n q 

n=l 3 

The three-state Potts model is defined by the partition function 



Z = exp [1 - Re(a*a i+(l )} (7) 

Wi} { «,A J 

where the variables Uj are defined on a three-dimensional hypercubic lattice 
and take the values 

o-i = exp ( — g- 1 J = 0,1,2 . (8) 

In the termodynamic limit the Potts model is known to undergo a (weak) 
first-order phase transition at /3 C = 0.36708(2) |L6| and the roughening tem- 
perature can be estimate to be (3 r ~ 0.93 (l^] . For (3 > f3 c the Z3 symmetry is 
spontaneously broken and the three ordered phases coexist, while at j3 = /3 C 
also the disordered phase coexists with the previous ones. 

In the finite cylindric geometry we are considering spontaneous symme- 
try breaking at low temperature cannot occur: the degeneracy of the ground 
state is removed, the energy of the symmetric, Z3 invariant, ground state 
being separated by an energy splitting E from the two degenerate mixed- 
symmetry states. 

The energy splitting is due to tunneling between the phases and is di- 
rectly linked to the free energy of the interface ||. According to the CWM, 
for r > (3 > p e , we assume 0, |, §, R > T, 



E(R, T) = C e~ aRT Z£ loop) ( |H • 4 2 loop) (R, T) (9) 



T 



C 



z (l loop) (1) 
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where 5 is an unpredicted constant and a convenient normalization has been 
chosen. 

We would like to stress that the two-loop contribution (|6|) does not 
depend only on the ratio z = R/T, like the one-loop term (g), but also on 
the minimal area A m = RT. If we put 2 = 1, A m = R 2 , in (||) we obtain 



E{R,R) =5 e~ aR2 jl+ ' 



where / = E2(i). The classical formula J 



£; d (i?)=Je-^ (11) 

can be recovered only neglecting the two-loop contribution. 

The comparison between formula (||) and the values of E extracted from 
MC simulations provides a simple and stringent way to verify the CWM 
predictions. 

It should be noted that no new free parameters are introduced within 
this approach: the formulae (|9|) and (|ll|) contain the same number of unde- 
termined parameters, namely a and 5. 

3 Monte Carlo results 

To extract the energy splitting E from MC-generated ensembles we follow 
the procedure of |]18|| . Defining the time-slice magnetization 

R T 

S'fc = — XI X] <r(zi,x 2 ,k) , (12) 

Xl=l X2 = l 

we compute the correlation function 

G(k) = (S S* k ) (13) 

where k = 0, 1, . . . , ^ , and we extract the transfer matrix low energy levels 
from the asymptotic fc-dependence of G(k) 

G(k)-Z = c {exp(-kE)+exp[-(L-k)E}} 

+c 1 {exp(-kE')+exp[-(L-k)E'}} + ... (14) 
Z = l + 2e~ LE + ... (15) 



(10) 
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Z = tr e being the partition function (the next-to-leading energy level 
E' turns out to be non-negligible in our range of parameters). 

Having so extracted the energy splitting E from the MC data for different 
values of the lattice sizes, we can compute the order-order interface tension 
a and the constant 5 by fitting our data with the formula (||). 

We have performed our simulations at /3 = 0.3680, the longest lattice size 
being fixed at L = 120, and the other sizes varying in the range 9 < T < 11, 
10 < R < 36, (R > T). This value of f3 is enough inside the ordered phase 
to make highly suppressed the probability of formation of order-disorder 
interfaces |ll|] but presents a correlation length large enough to make the 
lattice artifacts negligible and to consider domain walls as fluid interfaces. 

The fact that the disordered phase is substantially absent at this (3 can 
be seen from the histograms of the real part of the magnetization (l^] 

ReM = Re(^-J2 S ^j ( 16 ) 

as is shown, for example, in Fig.l. The modulus of the magnetization at 
this (3 is about 0.44 — 0.50 for the lattice sizes we are considering: in this 
figure the projection on the real axis of the ordered phases are clearly visible 
while the peak centered at ReM = 0, which would signal the presence of 
the disordered phase, is absent. 

Fig. 2 represents a typical distribution of the magnetization M for a 
sample of our MC-generated configurations. Most configurations consist of 
a single phase or of two phases separated by two interfaces (the minimum 
number compatible with periodic boundary conditions). The single-phase 
configurations are represented by the three clusters of points corresponding 
to the three degenerate vacua; the two-interface configurations form the 
straight lines joining these clusters. Three-interface configurations, which 
tend to fill uniformly the interior of the triangle, are clearly visible in Fig. 3: 
it corresponds to a T = 10, R = 20 lattice and to a larger probability of hav- 
ing tunneling events (i.e. interfaces), while these are much more suppressed 
for the lattice of Fig. 2 (T = R = 18). 

We have used a Swendsen-Wang cluster algorithm |2(| to perform our 
MC simulations. To keep under control correlations in MC time and cross- 
correlations between the G{k) observables, we have systematically scattered 
our measurements avoiding the measurement of two different observables at 
the same MC time. We have made between 0.6 • 10 6 and 1.8 • 10 6 sweeps for 
each experiment, depending on the lattice size, obtaining about 10 3 data per 
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observable. However, the covariance matrix turns out to be different from 
the diagonal form, which one expects from a sample of statistical indepen- 
dent data. We have taken this fact into account by including the covariance 
matrix in the fitting procedure to formula (|i~4| ) to extract the energy gap E: 
the results are reported in Tab. 1. The error on E has been estimated with 
an ordinary jackknife procedure. 
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Tab. 1. The values of E are reported with the x 2 P er degree of freedom 
and the confidence levels, as obtained from the fit of G(k) with formula 
(14)- The values in the last column are obtained from the best fit of all 
data to formula The same data are plotted in Fig. 4- 



7 



Fitting our results for the energy gaps E with the CWM formula (^) we 
obtain the following values of the interface tension and of the constant: 

a = 0.009912(75) 
5 = 0.1377(19) 

with a x 2 P er degree of freedom and a confidence level 

x 2 = 0.73 C.L. = 82% (17) 

thus confirming the accuracy of the CWM. In Tab. 1 the MC results for 
E are compared with the predictions of formula @ in which the best-fit 
values of 6 and a have been substituted. This comparison is represented 
graphically in Fig. 4. 

The importance of the inclusion of the two-loop contributions can be 
seen by fitting the MC data with the classical formula 

E(R,T) = 5 e~ uRT (18) 

and with the one-loop approximation {R > T) 

z (l loop) ( r\ 

E(R, T) = 5 e- aRT -^i ^ . (19) 

V > Z l q loop {\) V ' 

In the former case we obtain x 2 /d.o.f. = 36.3, in the latter x 2 /d.o.f = 3.60: 
the two-loop correction must be included to obtain a good agreement with 
numerical data. 

We have already noted that the two-loop corrections affect the value 
of E also for symmetric (T = R) lattices (cfr. (|l0|)), in contrast to what 
happens for the scale-invariant one-loop contribution (|5|) Q. Indeed, the 
importance of including two-loop corrections can be seen by fitting only the 
energy gaps E obtained on symmetric lattices (T = R): using the two-loop 
expression ( |To|) we obtain \ 2 = 0.79, while the classical formula ( |TT|) gives 
X 2 = 1.48. The results of all these fits are summarized in Tab. 2. 

We would like to stress the remarkable stability of the results obtained 
with the two-loop approximation fitting all the data or only the symmetric 
ones. This can be seen comparing the values of a and 5 given in the first 
and fourth line of Tab. 2. 
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0.01075(14) 


0.1866(36) 


1.48 


0.20 


T = 11 


2-loop 


0.00997(69) 


0.140(24) 


0.26 


0.90 


T = 11 


1-loop 


0.00965(69) 


0.151(26) 


0.22 


0.92 



Tab. 2. Results of the fit of E with two-loop, one-loop and classical 
approximations of the CWM, considering all values of(R,T), symmetric 
(R,R) lattices or (R,T=11) lattices. 



On the other hand, the result obtained using the classical formula (|lSj) is 
not compatible with the previous ones even using only the symmetric data, 
as it is shown in the fifth line of the same table. The best fit curve obtained 
from (18) in the latter case is plotted in Fig. 5 were also the "asymmetric" 
MC data are reported for comparison. 

We would also like to observe that a good agreement with the one-loop 
approximation of the CWM can be obtained Q if one considers low values 
of ratios z = R/T and high values of the minimal area A m = RT, i.e. where 
the two-loop contribution (JgJ) are maximally suppressed. This is show in 
the last two lines of Tab. 2 and in Fig. 6. 



4 Conclusions 

In this paper we have shown that the CWM in the two-loop approximation 
provides an excellent description of order-order interfaces in the 3D three- 
state Potts model. This result, together with the corresponding one for the 
3D Ising model [Q], strongly supports the hypothesis of the universality of 
the CWM description of interface physics in 3D statistical models. 

It is worth stressing again that the CWM corrections @ to the finite- 
size behavior of the interface free energy do not introduce any new free 
parameters with respect to the "classical" picture (|T§|). 

Besides the intrinsic physical interest of this picture, it should be noted 
that the CWM provides an accurate description of finite-size corrections 
to the free energy of rough interfaces, thus enabling one to extract correct 
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informations about physical observables from finite lattices of different ge- 
ometries, as it has been recently shown [21]. 
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Figure Captions 



Fig. 1. Histogram of the real part of the magnetization for a typical MC 
ensemble (in this case T = 20, R = 20). The absence of a peak in Re M = 
indicates that the disordered phase does not coexist with the ordered ones at 
our (3 = 0.3680. 

Fig. 2. Distribution of a sample of 5,000 configurations generated by the 
MC simulation in the complex plane of the magnetization for a lattice T = 
R = 18. The three clusters of points represent the one-phase configurations; 
the straight lines joining the clusters are the two-interface configurations. 

Fig. 3. The same of Fig. 2 for a lattice T = 10, R = 20. The 

three-interface configurations are uniformly distributed in the interior of the 
triangle. 

Fig. 4. Comparison of the predictions of the CWM with the MC data: 
the energy gap E is plotted as a function of z = R/T . The lines represent 
the best fit of all data to formula ^): from up to down they correspond 
respectively to T = 9, 10, 11 fixed (cfr. Tab. 1). 

Fig. 5. Comparison of the predictions of the classical formula with the 
MC data: E is plotted as a function of A m = RT. The line represent the best 
fit of the symmetric T = R data, reported with error bars, to formula jl^j); 
the asymmetric MC data are also reported: squares correspond to T = 9 
data, circles to T = 10 and diamonds to T = 11. 

Fig. 6. Comparison of the predictions of the two-loop (j^j and one-loop 
approximations for the T = 11 data, with z = R/T. The dashed line 
represent the one-loop best fit while the two-loop is the dotted line. 



12 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9310028vl 



This figure "figl-2.png" is available in "png" format from: 
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